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Circadian clocks exhibit the robustness of period and plasticity of phase against environmental 
changes such as temperature and nutrient conditions. Thus far, however, it is unclear how both are 
simultaneously achieved. By investigating distinct models of circadian clocks, we demonstrate reci¬ 
procity between robustness and plasticity: higher robustness in the period implies higher plasticity 
in the phase, where changes in period and in phase follow a linear relationship with a negative coef¬ 
ficient. The robustness of period is achieved by the adaptation on the limit cycle via a concentration 
change of a buffer molecule, whose temporal change leads to a phase shift following a shift of the 
limit-cycle orbit in phase space. Generality of reciprocity in clocks with the adaptation mechanism 
is confirmed with theoretical analysis of simple models, while biological significance is discussed. 

PACS numbers: 87.18.Yt, 05.45.Xt, 87.18.Vf 


Biological systems are both robust to external changes 
in the environment, and plastic to adapt to environmen¬ 
tal conditions. How are the robustness and plasticity, 
which seem to be opposing properties at a first glance, 
compatible with each other? In the present Letter, we 
address this question, by focusing on biological clocks, 
which are ubiquitous in organisms. 

Such biological clocks often work as pacemakers, to 
adapt to periodic events. One of the most prominent ex¬ 
amples of such oscillators is a circadian clock 0,0. To 
respond to periodic events, the following two criteria are 
generally imposed on a biochemical oscillator. 

1. Robustness of period: If the period of an oscilla¬ 
tor strongly depends on external conditions, the oscillator 
would not accurately predict time. For example, if the 
period of a circadian clock is sensitive to temperature, 
the clock malfunctions depending on the temperature. 
To avoid such error, the period of pacemakers should not 
be affected by external conditions such as temperature 
and nutrient compensation 0 0 . 

2. Plasticity of phase: The period of the circadian 
clock of most organisms is known not to correspond pre¬ 
cisely with 24 hours 0, and biological clocks are en¬ 
trained with the external 24-hr cycle [6| , so that the phase 
difference between the two does not increase with time. 
This entrainment is also necessary to adapt an abrupt 
change in the environment that may cause temporal mis¬ 
alignment between the internal and external cycles. For 
such entrainment, plasticity of the phase of the internal 
clock against external stimuli, e.g., changes in tempera¬ 
ture and/or brightness, is needed. 

Indeed, biological clocks satisfy both robustness and 
plasticity to changes in factors such as temperature and 
nutrient conditions, which change in the daily cwle. For 
examplejCircadian clocks of in vivo Drosophila [7|, Neu- 
rospora Q, and in vitro cyanobacteria 0, |l3] show tem¬ 
perature compensation of a period and are entrained by 
cyclic temperature changes. Robustness of period is also 
important to stable entrainment since it can reduce the 


difference between the period of inner clock and external 
cycle. In spite of some studies discussing the compatibil¬ 
ity between the two properties 0, however, little 

is known about the quantitative relationship between the 
two properties. 

To answer how the robustness of the period and 
plasticity of phase are compatible with each other, we 
first analyze two major models of a circadian clock, 
i.e., post-translational oscillator (PTO) SHEi and 
transcription-translation-based oscillator (TTO) [l2|, Il3l . 


15l| . which consists only of protein-protein interactions 


and both transcription and translation processes, respec¬ 
tively. Without imposing any special mechanism, we 
demonstrate that biological clocks with robustness of pe¬ 
riod against changes in an environmental factor generally 
exhibit phase entrainment against the cyclic change of 
that factor — reciprocity between the robustness of pe¬ 
riod and plasticity of phase: the plasticity increases with 
robustness. 


For PTO model, we adopt the KaiC allosteric model 
id, for in vitro cyanobacterial circadian clock system 
9j. Here, KaiC protein consists of six monomers, each 
of which has a phosphorylation site. The protein has 
active and inactive forms. Active (inactive) KaiC are 
phosphorylated (dephosphorylated) step by step, respec¬ 
tively. Phosphorylation reactions are facilitated with 
KaiA as an enzyme and dephosphorylation reactions 
spontaneously progress without an enzyme, kp and kdp 
denote the rate of phosphorylation and dephosphoryla¬ 
tion of KaiC, respectively, which depend on temperature 
as kp oc exp(—/JAp) and kdp oc exp(—fiEdp), where Ep 
{Edp) is the activation energy for phosphorylation (de- 
phosphorylatiion), respectively, and /3 is the inverse tem¬ 
perature by taking the Boltzmann constant as unity. The 
temporal evolution of the concentration of each phospho¬ 
rylated active (inactive) KaiC is given by rate equations 
(see model equations and Fig.lA of [3). 

This model shows a limit-cycle attractor in which the 
total phosphorylation level, i.e., the ratio of phosphory- 







2 



FIG. 1. Reciprocity between the robustness of period and 
plasticity of phase in the PTO model. (A) Difference between 
periods at two temperatures (/3i = 1.0 and /?2 = 1.5) {AT/T, 
red circle) and the amplitude of the phase response curve 
against a transient jump of temperature from /3i to /?2 {A<j), 
green square) plotted against different values of Edp while Ep 
is fixed at 1.0. AT is normalized by the period at /3i, and 
A(j) is normalized by the duration of stimulus and difference 
between /3i and 132 ■ ATjT and Ac)) are negatively correlated 
across the entire range of Edp- (B) Entrainability is plotted 
against various Edp with fixed Ep — 1.0. 

(C) Phase response curve against transient increase in /3. As 
a stimulus, the inverse temperature j3 is increased from jSi to 
/32 for the duration of one unit of time. Lines of different 
colors represent the PRCs for different values of activation 
energy for the dephosphorylation reaction Edp- 


lated monomers, oscillates in time. We demonstrated 
that the robustness of the period against various envi¬ 
ronmental changes is achieved by enzyme-limited com¬ 
petition [mH: With the increase in temperature, the 
abundance of the active form of the KaiC molecule in¬ 
creases, which in turn decreases the abundance of the 
free KaiA molecule, and thus the increase in the rate of 
phosphorylation kp is canceled out, when the total KaiA 
amount, Atotah is sufficiently small. This robustness in 
the period is achieved when Ep is sufficiently larger than 
Edp- We use the difference in periods between two dif¬ 
ferent temperature conditions (AT /T) as an indicator of 
the robustness of period. Its dependence upon Edp — Ep 
is given in FiglHA. 

This clock, on the other hand, entrains against exter¬ 
nal periodic change, so that the phase of phosphorylation 
oscillator coincides with that of external cycle. By im¬ 
posing external periodic change in temperature, we com¬ 
puted how many number of cycles are needed for the 
clock to entrain with this external cycle, and defined 
entrainability as the inverse of the number (see [ 13 ). 
Dependence of the entrainability and AT/T upon Edp 
with fixed Ep is plotted in FiglHA (red circle) and B. As 
Edp — Ep is smaller, AT/T becomes smaller and the en¬ 
trainability is higher. In other words, if the period of 


the clock is more robust against temperature change, it 
is entrained faster with the external temperature cycle, 
i.e., the phase has higher plasticity. 

Although this demonstrated the correlation between 
period robustness and phase plasticity, the entrainabil¬ 
ity here is a complicated indicator for the latter, as it 
can depend on the form of external cycle. Hence, we in¬ 
troduce a more tractable indicator for the plasticity of 
phase, by using a phase response curve (PRC) [2^. PRC 
is a function of phase and represents a phase shift intro¬ 
duced by a transient stimulus. When a transient stimulus 
is added to an oscillatory system, the period of oscilla¬ 
tion is temporally altered depending on the phase when 
the stimulus was added. The period finally returns to its 
original value. In this time, the phase of the oscillator 
progresses (or is delayed) from the original phase because 
of the temporal shortening (or lengthening) of the period. 
PRC represents such a phase shift Atp as a function of the 
phase (j) when the stimulus is applied. We computed PRC 
by transiently changing the inverse temperature from jSi 
to /32 for one time unit (see FiglT]C), by defining the 
phase of oscillation by the time when the total phos¬ 
phorylation level takes maximum at ^ = 0, 27r, • • •. As 
an indicator of the plasticity of phase, we measured the 
difference between maximum and minimum values of the 
phase change Acj) in PRC 2l| normalized by the magni¬ 
tude of a stimulus by fixing its duration as one time unit. 
The dependence of Acj) and AT/T on Edp with fixed Ep 
is plotted in Fig[T]4. When Edp is low, i.e., when the 
temperature dependence of dephosphorylation reaction 
is weak, AT/T is small and A0 is large. This reciprocity 
was also obtained against changes in other parameters, 
/3i and Atotai (see Fig.3 of (13). This indicates that a 
biochemical oscillator with a homeostatic period against 
an environmental change can easily shift its phase un¬ 
der the same environmental change. We also confirmed 
such reciprocity against change in ATP, i.e., the case of 
nutrient compensation (see Fig.5 of [13). 

Now, we examine if such reciprocity holds in the other 
class of circadian clocks, the TTO. In the TTO model, 
a clock-related gene is first transcribed and translated, 
and later such a translated protein represses the expres¬ 
sion of its own gene with a time-delay. When the tran¬ 
scription rate decreases, the amount of such protein also 
decreases, which weakens the suppression of the clock- 
related gene expression. Consequently, such genes are 
transcribed again, leading to the oscillation of the gene 
expression level. As a typical example of the TTO model, 
we choose here a model of a circadian clock of a fruit fly 
[^ (see model equations and Fig.IB of [ 13 ). By varying 
the activation energy for mRNA degradation, Ea, and 
fixing activation energies for other reactions, we mea¬ 
sured Acj) and AT/T using the same procedure as in the 
Kai model. Then, AT/T is low and A^ is high for a low 
Ea value, and AT/T {Acj)) increases (decreases) with the 
increase in Ea (FigUl and see also Fig.7 of [l7|). Thus, 
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FIG. 2. Reciprocity between the robustness of period and 
plasticity of phase in the TTO model. Difference between pe¬ 
riods at two temperatures (/3i = 0.0 and P 2 = 0.5) (AT /T, 
red circle) and the amplitude of the phase response curve 
against a transient jump of temperature from j3i to j32 (Ac)), 
green square) plotted against various Ea, while activation en¬ 
ergies for other reactions are fixed at 1.0. AT/T and A(j> are 
calculated similarly to how they are calculated in Fig[T] 



the reciprocity holds also in the TTO. 

To discuss the reciprocity analytically, we then study 
the Stuart “Landau model, a minimal model for simple si¬ 
nusoidal oscillation 23| . The model consists of the ampli¬ 
tude R and argument 0, where R and 0 reach a constant 
value at the limit-cycle attractor. Indeed, this model is 
derived as a normal form close to the Hopf bifurcation 
point. We introduce an external parameter /3: 


^^=/i(/?)i?-i?3, (la) 

^^ = /l(/3)a; + /2(^)i?^ (lb) 

where, /i(/3) is a response function of the first order term 
in complex Ginzburg-Landau equation, and /2(/3) is that 
of the third order term (for choice of each functions, see 
[l7|). Considering the stability of limit cycle, the re¬ 
laxation of R after perturbation is assumed to be much 
faster than that of 0 [ 2 ^. Here, the period is given as: 


T(/3) = 2^{/i(/3)(cu + /2(/3))}-'. (2) 


Thus, after the change (3^/3 + A/3, the dependence of 
period on (3 is given as: 

Alnr(/3) ~ -Aln/i(/3) - Ah{P) {uj + /2(/3))-'. (3) 


Here, we neglected higher order terms of A/3, assuming 
that it is sufficiently smaller than /3. From Eq. ©, 
if Aln/i(/3) = -A/2(/3 )(c. + /2(/3))-' (i.e., /{//i = 
~/ 2 /(^+/ 2 )) is satisfied, the dependence of the period on 
/2(/3) will be counterbalanced by /i(/3), and the period 
is compensated against a change in /3. 

The argument 0 is defined only on a limit-cycle orbit, 
and we introduce the phase </> to extend the definition to 
the phase space out of the limit-cycle attractor, in par¬ 
ticular to its basin. It is postulated that (p agrees with 
0 on the limit-cycle orbit, i.e., different orbits from the 
same cp converge to the same point on the limit cycle hav¬ 
ing the same 0 value. Now, we will derive an isochrone. 


which is a set of points with the same p on the phase 
space, p is expected to have rotational symmetry, hence 
the isochrone of the Stuart-Landau equation against the 
parameter /3 is derived as 

</)(i?,0,/3) =0 + /2(/3)|lni?-iln/i(/3)|. (4) 

(see a supplemental text of So Then, we consider an 
operation that increases /3 from /3o to /3o-l-A/3 and instan¬ 
taneously reverses it to /3o. By assuming that R instan¬ 
taneously relaxes to R*{l3o + A/3) = (/i(/3o -I- A/3))^/^ 
while 0 remains unchanged, the phase after the above 
operation is derived as 

<^(/3o+A/3) = ei(3o)+^^ {ln/i(/3o + A/3) - ln/i(/3o)} . 

(5) 

Hence, when A/3 ^ /3, the change in phase is derived as: 

A(/)(/3o) = /2(/3o)Aln/i(/3)/2. (6) 

Therefore, from Eqs. © and (IT6l) . changes in the period 
and phase are represented by an equality. 


a A In r -|- Ap = c, 


(7) 


where a = c = -/2(/3)A/2(/3)/2 (w-b/ 2 (^)), 

which depend only on /2(/3) and not on /i(/3). Thus, 
when we construct /i(/3), which compensates for the de¬ 
pendence of /2(/3) on /3 according to Eq. ([3|), the phase 
is altered as Ap = c. On the other hand, when /i(/3) is 
independent of /3, the phase is also independent due to 
Eq. (|T6l) . while the period is strongly dependent on (3 as 
AlnT = c/a = -A/2(/3)/(u; -b /2(/3)). 

We also confirmed the reciprocity is valid in the mod¬ 
ified van der Pol oscillator |25l | with strong nonlinearlity, 
i.e., beyond the neighborhood of Hopf bifurcation (See 
Eig.9 of 0). 

The origin of reciprocity is also understood from the 
viewpoint of adaptation motif. The standard minimal 
feedforward motif for adaptation consists of two compo¬ 
nents, X and y (26|. In the feedforward network in Figl3]4., 
an input changes both components x and y, while y gives 
an input to x. Here, the direct path to x and the indi¬ 
rect path via y from the input have opposite signs. Then, 
the response of the output x via the direct path is later 
canceled by y, and the adaptation behavior against the 
input is shaped. The degree of adaptation depends on 
the strength of the indirect regulation; weak regulation 
induces a partial adaptation and strong regulation lead¬ 
ing to the cancellation of the two paths, induces perfect 
adaptation [illii]. 

Our Stuart-Landau model also has a feedforward mo¬ 
tif consisting of amplitude and angular velocity. When 
an environmental condition /3 is changed, the angular 
velocity and amplitude are altered by the terms /i(/3) 
and /2(/5). After a direct change in angular velocity, the 
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FIG. 3. Schemes of the reciprocity between the robustness of 
period and plasticity of phase. (A) Schematic networks of a 
generic (bio)chemical oscillator exhibiting homeostasis of pe¬ 
riod. Pointed and flat arrowheads indicate positive and neg¬ 
ative regulation, respectively. Correspondences with a simple 
feedforward adaptation motif are represented by green charac¬ 
ters in parentheses. (B) Scheme of limit-cycle orbits with com¬ 
pensation of the period against environmental change. Blue 
dotted line is a stable limit-cycle orbit before environmental 
change. Green dashed line and magenta solid line are stable 
limit-cycle orbits after environmental change when the period 
is perfectly compensated and partially compensated, respec¬ 
tively. 


change is relaxed by the change in amplitude. The pe¬ 
riod is determined as the inverse of the angular velocity. 
If changes in the amplitude are large, period is perfectly 
compensated and phase is plastic. In contrast, if the 
change in amplitude is small, the angular velocity shows 
partial adaptation leading to partial compensation of the 
period while the phase is only slightly altered. Therefore, 
the reciprocity is understood as the adaptation dynamics 
on a limit cycle. 

Indeed, the above argument of the adaptation on the 
limit cycle generally holds, for PTO and TTO models, 
where we can generally consider the scheme of Figl3J4. 
Environmental change directly influences the angular ve¬ 
locity while it is also buffered in the amplitude and then 
influences the phase. In a biochemical clock, the period 
mainly depends on the rate-limit reactions, which are 
slower than others. Environmental change will alter the 
speed of such rate-limit reactions, which is later coun¬ 
terbalanced by the change in the concentration of buffer 
molecules [ 2 ^. In fact, in the PTO model, the amount 
of free enzyme working as a buffer molecule can counter¬ 
balance the speed of the rate-limit reaction. Hence, the 
period of the oscillator is homeostatic against environ¬ 
mental changes. Likewise, in TTO model, mRNA plays 
the role of such buffer molecule. 

In this time, the limit-cycle orbit of oscillators with 
compensation shifts in the phase space of chemical con¬ 
centrations to change that of a buffer molecule (see 
FigEB). When homeostatic response is achieved, the con¬ 
centration of a buffer molecule x should be changed with 
Aa; by the change in the external environment. Then 
the limit-cycle orbit will be shifted to change the con¬ 


centration of a buffer molecule, and the magnitude of 
such shift and the change in isocline will be 0(Aa:) con¬ 
sidering that continuous change in the isocline against 
Ax which is small. Then, A^ oc Ax is expected. On 
the other hand, when the change in the concentration 
of a buffer molecule is not sufficient to counterbalance 
the environmental stimulus, the concentrations of other 
molecules will change. Let us represent the concentra¬ 
tion of X needed for perfect adaptation as Ax*. Then, 
the change in the concentration of the other molecule of 
the lowest order is proportional to Ax* — Ax. The period 
also changes accordingly, so that AT/T oc Ax* — Ax is 
expected. By combining the two proportionally relation¬ 
ships, we obtain aA(j) + bATjT = Ax* with coefficients 
of proportionality a and b. 

We have shown that reciprocity exists in both the 
PTO and TTO models. The currently known mecha¬ 
nisms of circadian oscillation can be classified into the 
above two cases [l^, and the reciprocity is expected 
to be achieved universally in circadian clocks |30|. In 
a circadian clock system of a mold, Neurospora crassa, 
it was reported that a loss-of-temperature-compensation 
mutant, frq-7, shows smaller phase shift against transient 
temperature change than the wild type [8|, l3l|, l32|. Al¬ 


though the quantitative relationship between tempera¬ 
ture compensation and phase plasticity was not inves¬ 
tigated therein, we expect that a quantitative exper¬ 
iment will confirm our reciprocity, not only in Neu¬ 
rospora crassa but also in other organisms in which loss- 
of-temperature-compensation mutants are isolated, e.g., 
fruit fly [33| and cyanobacteria [3^. Here, we demon¬ 
strated the reciprocity against changes in the tempera¬ 
ture and the nutrient concentration, but from theoreti¬ 
cal consideration, it is expected to hold generally against 
a variety of stimuli, such as the change in strength of 
light and transcription rate [illai, as long as the adap¬ 
tation mechanism works. Moreover, it is also expected 
that the reciprocity is not limited to the circadian clock; 
it holds generally as long as the adaptation mechanism 
with buffering molecules works 37| . Our reciprocity will 
give a general quantitative law for such adaptation sys¬ 
tems. 
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SUPPLEMENTAL MATERIAL 

MODELS 


Post-translational oscillator (PTO) model 


We introduce the KaiC allosteric model [1^ (Fig. 
UK). The KaiC protein has six monomers, and each 
monomer has multiple phosphorylation sites. Here, we 
assumed each KaiC monomers have only two phospho¬ 
rylation states, phosphorylated and unphosphorylated. 
KaiC hexamer takes an active or inactive form. By de¬ 
noting Ci and Ci as active and inactive forms with the 
i phosphorylated monomers, respectively, their temporal 
changes are given as: 


d[Ci\ 

dt 


d[Ci] 

dt 


(1 — 


kp[A][C,.^] 
Ki—i + [A] 


(1 - difi) 


kpjAm] 

Ki + [A] 


+S^,ob[Ci]-S,,ef[a], 


(la) 


kdp{{l - S,,e)[a+i] - {1 - S,,o)[a]) 

-S,,ob[Ci]+S,,ef[a], (lb) 


A 


total — 






[^][a] 

Ki + [^] 


(Ic) 


where A denotes the free KaiA protein that works as 
an enzyme for phosphorylation. Atotai is the total KaiA 
amount, which is a constant because the total amounts 
of both KaiC and KaiA are conserved quantities, and 
[a;] denotes the concentration of a;. ATj is the dissociation 
constant between Ci and A. kp and k^p denote the rate of 
phosphorylation and dephosphorylation of KaiC, respec¬ 
tively, which depend on temperature as kp oc eypp{—j3Ep) 
and kdp oc exp{—j3Edp), where Ep {Edp) is the activation 
energy for phosphorylation (dephosphorylatiion) and P 
is the inverse temperature by taking the Boltzmann con¬ 
stant as unity. For case of the nutrient compensation 
(Fig. [S]), we considered the model where only the phos¬ 
phorylation reaction speed, kp^ is proportional to ATP- 
to-ADP ratio and others are independent of it. 


Transcription-translation-based oscillator (TTO) 
model 


We introduce a model of a circadian clock of a fruit 
fly 01 (FigHK). The governing differential equations are 


A AC, AC, AC, AC, AC, AC, 

tl\tl\tl\tl\ll\ti\ 



It 


(j) < - M -► R 


FIG. 1. Schemes of models of the circadian clock. (A) KaiC 
allosteric model as an example of the PTO model. Ci and 
Ci are an active and inactive substrate (KaiC) with i phos¬ 
phorylated residues, respectively. A is an enzyme (KaiA) and 
ACi is an enzyme-substrate complex with i phosphorylated 
residues. (B) TTO model. M is an mRNA and R is the 
precursor of a protein. Q and P are an extranuclear and 
internuclear protein, respectively. P can negatively regulate 
gene expression and create a negative feedback loop. 


described below: 


d[M] 

k 

a[M] 


(2a) 

dt 

h+[P] 

a' + [M] ’ 


d[R] 

s[M] 

b[R] 

c[Q] 

(2b) 

dt 

s' + [M] 

b' + [P] 

d -b [Q] 

d[Q] 

b[R] 

c[Q] 

d[Q] 

(2c) 

dt 

V + [R] 

d -b [Q] 

d' + [Q] 


u[Q] v[P] 
u' + [Q] v' + [P] ’ 

d[P] _ u[Q\ v[P\ 

dt u' + [Q] v' + [P] ’ 


where M is the mRNA of the clock-related gene {per 
mRNA); R is the protein precursor of a clock-related pro¬ 
tein (PER protein); Q and P are an extranuclear protein 
and a nucleic protein, respectively; and [x] denotes the 
concentration of x. fc, a, s, b, c, d, u, and v are rate con¬ 
stants, and h, a', s', b', c', d', u', and v' are dissociation 
constants. 


In 


22| , it was reported that a and k are especially im¬ 


portant to determine the length of the period. Following 
this report, we set the rate constant of each reaction to 
follow the Arrhenius equation, i.e., a kinetic constant of 
mRNA degradation as a oc exp(—/3Po) and that of tran¬ 
scription as fc oc exp(—/3Pfc) where Ea and Ek are acti¬ 
vation energies of mRNA degradation and transcription, 
respectively. 
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60 


FIG. 2. Entrainment of the clock in the PTO model against 
external temperature cycle, For different values of activation 
energy of dephosphorylation reactions Edp- (A) Time evo¬ 
lution of phosphorylation level {Ei[Ci\/f)Ctotal) with Edp = 
0.1 (cyan line), 0.4 (green line), and 0.6 (orange line), plotted 
against the time, normalized by the period of temperature 
cycles. As shown with the red line, the temperature (to be 
precise, its inverse /3) is periodically changed between — 1.0 
and = 1.1 with the period at fii, where each interval at 
/3i and /I 2 is set identical. The phosphorylation oscillation 
is entrained with the external temperature cycle, at the time 
shaded in the figure. (B) Plots of phosphorylation levels per 
period of the external cycle, i.e., at the time when the temper¬ 
ature is switched from /3i to ^ 2 - Initially, the phosphorylation 
level changes per period, and then, after entrainment, phos¬ 
phorylation level takes an almost constant, when the clock 
is entrained with the external cycle. For smaller Edp, the 
phosphorylation clock is entrained faster. 


CALCULATION OF THE ENTRAINABILITY 


To calculate the entrainability in Fig.IB in the main 
text, initially, 36 oscillators are set at same intervals of 
the phase, and the number of cycles needed for the oscil¬ 
lators to synchronize is computed. As for the numerical 
criteria, and we regard that the clock is entrained if the 


A 0.45 


0.0 

■0.05 
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> 
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FIG. 3. Reciprocity between the robustness of period and 
plasticity of phase in the PTO model under various enzyme 
concentrations and under various temperatures. (A) Differ¬ 
ence between periods at two temperatures (/3i = 1.0 and 
/32 = 1.5) {AT/T, red circle) and the amplitude of the phase 
response curve against a transient jump of temperature from 
Pi to P2 {Ap, green square) plotted against various total con¬ 
centrations of the enzyme. (B) Difference between periods 
at two temperatures (/3i and P2 = Pi + 0.5) (AT/T, red cir¬ 
cle) and the amplitude of the phase response curve against 
a transient jump of temperature from Pi to P 2 {A(j>, green 
square) plotted against various Pi. ATjT and Ap are calcu¬ 
lated similarly to how they are calculated in Fig.l in the main 
text. 


circular variance of phase from different initial conditions 
is smaller than 10“®. Here, the temperature cycle is ap¬ 
plied as a square wave between Pi = 1.0 and P 2 = 1.1 
with an equal interval, with the period at the condition 
of P = Pi- When entrainability is zero, the oscillator is 
never entrained. See also Fig.2. 


ANALYSIS OF STUART-LANDAU EQUATION 


We introduce the Stuart-Landau equation with an ex¬ 
ternal parameter /3: 

=/i(/3)i?-i?3, (3^) 

at 

^^=f,(P)uj + f2iP)R^. (3b) 

This form is derived from the complex Ginzburg-Landau 
equation, where /i(/3) represents the change in the bifur¬ 
cation parameter, and f 2 {P) in that for phase-amplitude 
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FIG. 4. Phase response curve of the PTO model against the 
transient increase in /3 for different values of (A) Atotai and 
(B) /3i, the inverse of temperature. PRCs are calculated in 
the same manner as in Fig.lC in the main text. (A) Lines of 
different colors represent PRCs for different values of the total 
KaiA Atotai', Atotai = 0.01 (red line), 0.015 (orange line), 0.02 
(yellow line), 0.025 (green line), 0.03 (cyan line), 0.035 (indigo 
line), and 0.04 fiM (purple line). (B) Lines of different colors 
represent PRCs for different values of /3i, while the inverse 
temperature is changed transiently to /32 = /3i + 0.5; /3i = 0.0 
(red line), 0.1 (orange line), 0.2 (yellow line), 0.3 (green line), 
0.4 (cyan line), and 0.5 (purple line). 


coupling. 

HA 

— = him + ^lv)A - (1 - ^MP))\A\^A, (4) 


0-01 ( gM ) 0.03 

FIG. 5. Reciprocity between the robustness of period and 
plasticity of phase in the PTO model against the change in 
ATP ratio ([ATP]/([ATP] + [ADP])) and enzyme concentra¬ 
tions. (A) Difference between periods {AT/T, red circle) for 
two ATP ratios (initial ATP ratio x and a; -I- 0.1), and the am¬ 
plitude of the phase response curve against a transient jump 
of ATP ratio from a; to a: -|- 0.1 (A/), green square) plotted 
against different values of initial ATP ratio x. (B) Difference 
between periods {AT/T, red circle) for two ATP ratios (x 
and a: -I- 0.1), and the amplitude of the phase response curve 
against a transient jump of ATP ratio from a; to a: -f 0.1 {A(/>, 
green square) plotted against different values of total concen¬ 
trations of the enzyme. Here, we set initial ATP ratio x as 
0.5. AT /T and A(/> are calculated similarly to how they are 
calculated in Fig.l in the main text. 


where, the first order term (/i(/3)) and the third order 
term (/2(/3)) should have different dependency upon /3 
for the adaptation mechanism to work. Indeed the reci¬ 
procity holds generally for other forms of the third order 
term to satisfy adaptation. 

Considering the stability of limit cycle, the relaxation 
of i? after perturbation is assumed to be much faster than 
that of 0. Hence, R*, the steady-state value of R (> 0), 
is obtained as 


R*i/3) = (/i(/3))'/^ 

(5) 

Here, the period is given as: 


T(/?) = 2^{/i(/3) {oj+ hmr^ ■ 

(6) 

Here, the argument 0 is defined only on a limit-cycle 
orbit, and we introduce the phase (j) to extend the defini- 


tion to the phase space out of the limit-cycle attractor, in 
particular to its basin. It is postulated that </> agrees with 


0 on the limit-cycle orbit, i.e., different orbits from the 
same (/> converge to the same point on the limit cycle hav¬ 
ing the same 0 value. Now, we will derive an isochrone, 
which is a set of points with the same (j) on the phase 
space. (/) is expected to have rotational symmetry and is 
given as 

^{R,e,l3) = e{/3)+giR{/3)). (7) 

Moreover, the time evolution of (j) should coincide with 
that of 0. Thus, 

^ = f,{/3)ut +hi/3) M/3). (8) 

From Eqs. (EiD, (|3b|, (0), and (IS]), the time evolution of 
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FIG. 6. Phase response curve of transcription-translational- 
based oscillator (TTO) model against the transient increase in 
P for different values of if a- The zero-phase point 0 = 0, 27r is 
defined as the state in which the total amount of protein ([i?] -|- 
[Q] + [-P]) takes its maximal value, and the phase increases 
from 0 to 27r proportionally to time. The inverse temperature 
P is increased from Pi = 0.0 to P 2 = Pi + 0.5 = 0.5 for the 
duration of one unit of time. Lines of different colors represent 
PRCs for different activation energies for mRNA degradation 
Ea- Ea = 0.0 (red line), 0.2 (orange line), 0.4 (yellow line), 
0.6 (green line), 0.8 (cyan line), and 1.0 (purple line). 

0.3 
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1.0 F 2.0 

■-k 

FIG. 7. Reciprocity between the robustness of period and 
plasticity of phase in the TTO model. Difference between pe¬ 
riods at two temperatures {Pi = 0.0 and P 2 — 0.5) {AT/T, 
red circle) and the amplitude of the phase response curve 
against a transient jump of temperature from /3i to P 2 {Ap, 
green square) plotted against various Ek where activation en¬ 
ergies for other reactions are fixed at 1.0. ATjT and Ap are 
calculated similarly to how they are calculated in Fig.l in the 
main text. 



p is derived as: 

dp{R,&,P) dQ dR dg{R) 

dt dt dt dR 

= fi{P)uj + h{P)R^ 

+{/i(/3)i?-i?n^. 

Hence, 

{fi{P)R - = /2(/3){/i(/3) - R^}. 

Thus, 

dgjR) f2W) 
dR R 


(9) 

( 10 ) 

( 11 ) 


Accordingly, we obtain 

g{R) = f 2 {P)\nR + C. (12) 

Next, we will derive C. Because the phase p{R,Q,P) 
should coincide with 0 at i? = R*, g{R*) becomes zero, 
and then from Eq.(4) in the main text and Ea. dT^ . 

C = -i/2(^)ln/i(/3), 

0(i?,0,/3) = 0 + /2(^)|lni?-iln/i(/3)|. (13) 

Thus, the isochrone of the Stuart-Landau equation 
against the parameter /3 is derived. Then, we consider 
an operation that increases P from Po to /?o + A/3 and 
instantaneously reverses it to Po- At this time, by as¬ 
suming that R instantaneously relaxes to R*{Po + A/3) = 
ifiiPo + A, 5 ))^A while 0 remains unchanged, the phase 
after the above operation is derived as 

P{Po+Ap) = e{Po)+^^ {lii/i(/3o + A/3) - ln/i(/3o)} . 

(14) 

In contrast, the phase before the operation is given as 
P{Po) = 0(/3o). (15) 

Hence, when Ap P, the change in phase is derived as: 

A(^(^o) = <)>(/3o + A/3) - </.(^o) 

= f2iPo) {lii/i(/3o + A/3) - ln/i(^o)} /2 
= /2(/3o)Aln/i(/3)/2. (16) 

Therefore, from Eq.(6) in the main text and Ea. dlBL 
changes in the period and phase are represented by an 
equality. 

aA In T-f A^ = c, (17) 

where a = /2(/3)/2, c = -/2(/3)A/2(/3)/2 (w-b/2(/3)). 
The reciprocity is generally true also for some different 
forms, such as: 

= h{P)R-MP)R^ (18a) 

= /l(^)cc + /2(/3)i^^ (18b) 

= h{P)R-h{P)R^ (19a) 

= h{P)u + R^. (19b) 

In both the cases, Eq. in} still holds, while the expres¬ 
sion of a and c are different. Eor Eq.(18), a = 1/2 
and c = —Aln/2(/3)/2. For Eq.(19), a = l/2/2(/3) and 
c = {A/2(/3)(/2(/3)a; + l)-i - Aln/2(/3)}/2/2(/3). 


dR{P) 

dt 

dQ{P) 

dt 


dR{P) 

dt 

de{P) 

dt 
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FIG. 8. Phase response curve of the van der Pol oscillator 
against the transient increase in /3 for different values of Ei. 
The zero-phase point 0 = 0,27r is defined as the state in 
which X takes its maximal value, and the phase increases from 
0 to 27r proportionally to time. The inverse temperature /3 
increases from /3i = 0.0 to j32 = /li+0.5 = 0.5 for the duration 
of one unit of time. The strength of nonlinearity e is e = 0.1 
for (A) and e = 2.0 for (B). Lines of different colors represent 
PRCs for different values of E\\ E\ — 0.0 (red line), 0.2 
(orange line), 0.4 (yellow line), 0.6 (green line), 0.8 (cyan 
line), and 1.0 (purple line). 

VAN DER POL OSCILLATOR MODEL WITH 
PARAMETER 13 


FIG. 9. Reciprocity between the robustness of period and 
plasticity of phase in the van der Pol model. Difference be¬ 
tween periods at two temperatures (/3i = 0.0 and /32 = 0.5) 
(AT/T, red circle) and the amplitude of the phase response 
curve against a transient jump of temperature from to /I 2 
(A<(), green square) are plotted against various Ei with fixed 
E 2 = 1.0 in the van der Pol model with (A) e = 0.1 and (B) 
e = 2.0. ATjT and A(j!> are calculated in a similar way to 
Fig.l in the main text. 

where /3 is an environmental parameter. If e is sufficiently 
small, the amplitude A can be derived by using pertur¬ 
bation calculation as 


The van der Pol oscillator is one of the simplest 
nonlinear-oscillator models, and it is given by: 




( 20 ) 


The above equation can be decomposed into two ordinary 
differential equations as follows: 

dx 
dt 
dy 


= y, 


dt 


= e(l — x^)y — bx, 


(21a) 

(21b) 


where e is the strength of nonlinearity. As e increases, 
the system deviates more from a harmonic oscillator. 

Here, we modify the van der Pol oscillator to show the 
change in the amplitude against a change in an external 
parameter as in the Stuart-Landau equation. We alter 
van der Pol oscillator as 

dx . 

= y, (22a) 


dt 

dt 


= e(/i(/3)-a:^)y-6a;. 


(22b) 


dt V 2 8 y ’ 

(i^rf = /i(/3), 


(23) 

(24) 


where |A|* is a fixed point value of |A|. On the other 
hand, we introduce the dependence of the velocity on 
the environmental parameter as 


1 dx 
1 dy 


= y, 

= e(/i(/?) 


x^)y — bx. 


(25a) 

(25b) 


When e is small, the above modified van der Pol oscillator 
is expected to demonstrate same behavior as the Stuart- 
Landau model, as described in the main text. When e 
is large, however, the nonlinearity becomes large and the 
oscillatory behavior is altered from sinusoidal to relax¬ 
ation. 

To simulate the above model, we choose /i(/3) and 
/2(/3) as exponential forms similar to the Arrhenius equa¬ 
tion in biochemical oscillators, i.e., /i(/3) = and 
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/2/3 = e Then, the above equations are given as 


— = e ^ ^y, 
dt 

dt 


(26a) 

(26b) 


where Ei is given by = A + E 2 . We use the above 
equations. 

When the intensity of nonlinearity, e, is small, the 
magnitude of change in the period and magnitude of 


change in the phase are fitted well by a linear relation¬ 
ship aAT/T + bA(j) = c (a, b, and c are constants) (see 
FiginK). Here, as the intensity of nonlinearity increases, 
the orbit of a limit cycle is deformed from a circle, and 
the dynamics shifts from sinusoidal oscillation to relax¬ 
ation oscillation. Still, the reciprocity is valid as long 
as the magnitude of stimuli is not exceedingly large (see 
FiginiB). Thus, the reciprocity is a universal feature be¬ 
yond the neighborhood of Hopf bifurcation. 


